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We define the new description of the electromagnetic current to hold the current conserva- 
tion in the momentum-dependent Dirac fields from the Ward Takahashi identity. To describe 
the momentum-dependence we solve the relativistic Hartree-Fock approximation by using 
the one-pion exchange. In addition we discuss on contribution from the one-pion exchange 
current and the core polarization. It is shown that this current can reduce the convection 
current in the isovector case, whose value has been too big due to the small effective mass 
in the usual relativistic Hartree approximation. 
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I. INTRODUCTION 



The past decades have seen many successes in the relativistic treatment of the nuclear many-body 
problem. The relativistic framework has big advantages in several aspects 111: a useful Dirac 

HQ 

phenomenology for the description of nucleon-nucleus scattering the natural incorporation 

of the spin-orbit force Q] and the saturation properties of nuclear matter in the microscopic 
treatment with the Dirac Brueckner Hartree-Fock (DBHF) approach 

These results conclude that there are large attractive scalar and repulsive vector-fields, and 
that the nucleon effective mass is very small in the medium. However this small effective mass 
leads to small Fermi velocity, which makes some troubles in the nuclear properties: too big 
magnetic-moment [a] and too big excitation energy of the isoscalar giant quadrupole resonance 
(ISGQR) state 61. As for the isoscalar magnetic- moment, this enhancement is cancelled by the 

Q Q 

Ring-Diagram contribution p|; this relation is completely realized by the gauge invariance [g|. 
As for the isovector one, however, this contribution does not plays a significant role because the 
symmetry force is not efficiently large. 

In this subject most of people believed that the momentum dependence of the Dirac fields is 
negligible in the low energy region, particularly below the Fermi level. A momentum dependence 
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of the Schrodinger equivalent potential automatically emerges as a consequence of the Lorentz 
transformation properties of the vector-fields without any explicit momentum dependence of the 
scalar and vector fields. In fact, only very small momentum dependence has appeared in the 
relativistic Hartree-Fock (RHF) calculation P. llo|. 

In the high energy region, however, the vector-fields must become very small to explain the 
optical potential of the proton-nucleus elastic scattering P, [n|], and the transverse flow in the 
heavy- ion collisions [l^]- The momentum-dependent part is not actually small though it has 
not been clearly seen in the low energy phenomena. Furthermore S. Typel introduces the 
non-local parts and succeeded to improved nuclear properties. 

In the previous paper 1 14] we showed that the momentum dependence of the Dirac-fields is very 
sensitive to the Fermi velocity though it affects little the nuclear equation of state. In that work 
we introduced the one-pion exchange force, which produce the dominant contribution of the 
momentum dependence and suppresses the Fermi velocity, and explain the ISGQR energy. 
We can easily imagine that the one-pion exchange force largely produce the momentum depen- 
dence because the interaction range is largest. Since the momentum-dependent fields break the 
current conservation, we have to define the new current caused by the vertex correction. 
In this paper, thus, we investigate the nuclear current using the momentum-dependent Dirac 
fields. For this purpose we define a new current to hold the current conservation in the 
momentum-dependent Dirac fields, and discuss its effect on the the nuclear static current. In 
this work we focus only on the convection current which is sensitive to the Fermi velocity, and 
omit the spin current. 

In the next section we explain our formalism to make a conserved current under the momentum- 
dependent selfenergies. In Sec. 3 we show our numerical results for the static current in our 
formulation. Then we summarize our work in Sec. 4. 



II. FORMALISM 
A. Nucleon Propagator 

Now we describe the propagator of nucleon with momentum p in the isospin space as follows. 

sw-f™ ° ) (i) 

V S n (p) J 

Here we assume the spin isospin-saturated nuclear matter, and define the proton and neutron 
propagator as 

Sn(p) = S P ( P ) = Snip). (2) 
The nucleon propagator in the selfenergy £ is given by 

S N 1 (p) = i>-M-i;(p), (3) 
where has a Lorentz scalar part U s and a Lorentz vector part J7«(p) as 

Z(p) = -U s (p) + ^U fl (p). (4) 
For the future convenience we define the effective mass and the kinetic momentum as 

M*(p) = M-U s {p), 

n M (p) = p„ - u u ( P ). (5) 

The single particle energy with momentum p is obtained as 



e(p) = Po 



on— mass— shell 



= y / U 2 (p) + M* 2 + U (p). (6) 
Then the detailed form of the nucleon propagator eq. (jSJ) is represented by 

Sn(p) = S F (p) + S D (p) (7) 

with 

Sf(p) = {flip) + M *(P)} n 2_M*2 + t§ (8) 

= 2i7r{ nip) + M*(p)}n(p)8( Po )5[V(p)}, (9) 

where n(p) is the momentum distribution, and 

y(p)^i{n 2 (p)-M* 2 ( P )}. (io) 
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B. Momentum-Dependent Selfenergies 

We can easily suppose that it is the one-pion exchange force which produces the major momen- 
tum dependence because the interaction range is largest. In this work, thus, we introduce the 
momentum dependence to the Dirac fields due to the one-pion exchange, and discuss how the 
Fock parts affects the nuclear current. 

Along this line we define a Lagrangian density in the system as 



C = - M)if> + X -d^ a d^ a - \ml4> a $a ~ U[a] + \ml^^ 



f — - - C IV 

PitaPa + i^-^lbl^Ta^d^a + g^ipipa - g^l^^ - -^pi^l^} 2 (11) 
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where tp, (ft, a and u are the nucleon, pion, sigma-meson and omega-meson fields, respectively, 
and the suffix a indicates the isospin component. In the above expression we use the pseudo- 
vector coupling form as an interaction between nucleon and pion. The selfenergy potential of 



the (T-field U[a] is given as Ref. 
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3- 



U W] = \mW + \B^ + \C^ 

The symbols m n , m a and are the masses of 7T-, a- and w-mesons, respectively. In addition, 
we also introduce the isovector nucleon- nucleon interaction into the Lagrangian (jllj) . so as to 
discuss on the isovector current later. 

Next we calculate the nucleon selfenergies. The nucleon selfenergies are separated into the local 
part and the momentum-dependent part as U Q (p) = + U^(p), where a = s,/i. The a- and 
o;-meson exchange parts produce only very small momentum dependence of nucleon selfenergies 
0, ^| as their masses are large. In fact the RH and RHF approximations do not give any 
different results in nuclear matter properties after fitting parameters of a- and w-exchanges 
On the other hand the one-pion exchange force is a long range one and makes for a large 
momentum dependence while it does not contribute to the local part in the spin-saturated 
system. Subsequently we make the local part by RH of the a- and u;-meson exchanges, and the 
momentum-dependent part by RHF of the pion exchange, and thus we omit the kinetic energy 
part mesons except pion in ea- Hllf) . This method is shown in Ref. 
within the RHF framework. 



Ill ] to keep the self-consistency 
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In this model the local part of the selfenergies are given as 

= g*(<7) (13) 
Ui = do^PH (14) 

where (a) is the scalar mean-field obtained as 

^U[{a)]=g aPs (15) 

In the above equations the scalar density p s and the vector Hartree density ph are given by 

4 f d 3 p , ,M*(p) . . 

* = 4 i W? n{p) W (16> 

« " 4 f IW>BTT- < 17 > 
7 {2-kY n (p) 

where n(p) is the momentum-distribution, and ILj (p) is defined by 

n>) = ^[n 2 (p)-M* 2 fr)] (is) 

As a next step we define the momentum-dependent parts of the selfenergies as the Fock parts 
with the one-pion exchange. When using the pseudo vector (PV) coupling the Fock parts do not 
become zero at the infinite limit of the momentum \p\. One usually erases these contributions by 
introducing the cut-off parameter. In this work, instead of that, we subtract these contributions 
from the momentum-dependent parts (these contributions can be renormalized into the Hartree 
parts): U a — > U a — U a (p — > oo). Thus we obtain the momentum-dependent parts of the 
selfenergies as 

tF,^ 3/2 /• d 3 k _, t ,M*(fc) 



2 J (2vr)^ n (fc) 



T F, s . 3/2 f d 3 k f ,Mk) 



UZ( P ) = -^J^n(k) T ^A 7r ( P -k), (20) 
where the A^g) is the pion propagator defined as 

M?) = 2 1 2 - ( 21 ) 

In the above vector selfenergies we omit the tensor-coupling part involving [![(&)• (p— k)](p — k)^. 
This term is very small if the selfenergy is independent of momentum 9] , and their momentum 
dependence is actually very small as shown later. 
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C. One-Body Current Operator 

If the selfenergy has a momentum dependence, the current operator must be also changed to 
satisfy the current conservation. We then define the current vertex T fl (p + q,p) as 

r^p + q,p)=^ + A^(p + q,p). (22) 

The Ward-Takahashi (WT) identity gives the following relation about the current vertex 

S N (P + q)q^S N (p) = -S N (p + q) + S N (p). (23) 

This expression is rewritten as 

q"T^ = S] v 1 (p + q)-S^(j ) ). (24) 

Substituting eq.(jSJ) into the above equation (|24[). the density-dependent vertex correction is 
obtained as 

q^(p + q,p) = -Z(p + q)+Z(p). (25) 

In this work we restrict our discussion on the convection current, not on the spin current. In 
Appendix A we show an approximate method to derive the vertex correction from the above WT 
identity though the anomalous current, which is proportional to cr fJiU q 1 ', still has an ambiguity 
because the WT identity cannot make any condition for it. 

When we use the one-pion exchange force, we have to take into account the meson-exchange 
current. In Appendix B we show that our approximate formulation is also satisfied in the electro- 
magnetic current by including the one-pion exchange current, and the electromagnetic current 
operator is proportional to (l + T3)/2, namely only protons contribute to the convection current. 
In this work we focus on the static current in the nuclear matter, then we need to get the zero 
limit of the momentum transfer q. In this limit the vertex correction becomes 

A» = limA^ + ? ,p) = — E(p). (26) 
g^o dp^ 

Using the above vertex correction, the current density of the whole system is given as 

n, 



f d 4 p 1 + r 3 <9£ 



A r n«(p)^(p) - n»^M + M*(p) d -^l } S[V(p)] 



(2vr) 3 M dpa dp^ 



d3 P „&0(p) 5> } 



(2tt) 3 ^' n ( p ) 



(27) 

po=e(p) 
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where is the Fermi distribution function for proton, and is defined by 



= ^{n 2 W- m * 2 (p)}- (28) 

Let us consider the one-particle state on the Fermi surface. The space current density contributed 
from this nucleon can be written as 

3 = f^j\\\p\=p F = D p e(p)\ lpl=PF , (29) 

where the total derivative D p is defined on the on-mass-shell condition : pp = s(p)- The above 
equation completely agrees with that derived by the semi-classical way 
In the non-relativistic framework the effective mass is defined by 

A^ = (2^e(p)r 1 ||p|- P „ (30) 
which is so called the 'Landau mass'. Then the above spatial current density is 

3 = ^. 



(31) 



In our case including the momentum-dependent Dirac fields, the value of the Landau mass M£ 
cannot be uniquely determined from the relativistic effective mass M* while in the Hartree 
approximation the Landau mass becomes M£ = IIo(pf)- 



D. Core Polarization Current 

As for the actual nuclear current observed by experiments the core polarization also plays an 
important role. Of course this effect cancels contribution of the effective mass in the isoscalar 
case [it] ]. 

Here we should consider a system that one valence nucleon populates a state on the Fermi surface 
of the saturated nuclear matter. In this system the momentum-distribution can be described as 

n{p, t) = no(p) + ^An(p, r), (32) 

where no(p) = 6(pf — \p\) shows the usual Fermi distribution with the Fermi momentum pp, 
and An(p) cx S(\p\ — pp) indicates the valence nucleon part. The suffix r decades the isospin 
for the valence nucleon, 
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The valence nucleon varies the selfenergies of nucleons below Fermi surface from that at the 
saturated matter as 

U a (p)^U a (p) + AU a (p). (33) 
In addition the function V(p) is also varied as 

V(p) = V (p) + AV(p) (34) 

with 

AV(p) = -IPA£/ M + M*AU S . (35) 

The current density is described with the following expression. 

■n ST f d *P ft \ dV (p) f oa\ 

^ = -^J (2^ /(p ' T) ^' (36) 

where f(p, r) is the four-dimensional momentum-distribution for nucleon with isospin r, which 
is given as 

f(p,r) = n(p,T)S(p - e p ) (37) 
= J—n(p,T)6[V(p)}e(p ). (38) 

n (p) 

The variation along eqs. ()32|) — ()35|) leads to the above four dimensional momentum-distribution 

f(p,r) as 

f(p,r)=Mp,T) + Af(p,r) (39) 

with 

fo(p,r) = n (p)5[V (p)]e(p ), (40) 
Af(p,r) = An(p, r)8[V o (p)}0(p o ) + n (p) | j AF(p,<r)0(po). (41) 

The first term comes from the valence nucleon, and the second one from the core polarization. 
The total current density is given as 



^PB+f val +fcor- (43) 



& = / l^fmVip)]^ (42) 
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The first term is the current density of the saturated matter, and the second current density 



shows the contribution from the valence nucleon as 



J val 



/^ An(p)i(w _ Sp) |M. m 

The third current density jcor is so called the core polarization current density, which is caused 
by the variation of the selfenergies of nucleons in Fermi sea and given by 

fcor = -2^/ T^MP^)[^ r m{p)\ + -^ r {^^}v=vAV(p)] (45) 
= " 2 E / (03"o(p,r)A{AF(p)5[y o (p)]}. (46) 

Here it should be noted that the time component of the core polarization current density does 
not change the nucleon density: 

fcor = " 2 E / -03 n o( P ,r)-^{AV(p)S[V o (p)}} = 0. (47) 

Now we define the z-axis as the direction of the current at the matter. First we calculate the 
isoscalar current density by taking the valence nucleon part of the momentum distribution to 
be 

An(p,T) = ^f5(p-a) (48) 

with 

a = p F z, (49) 

where Q is the volume of the system, which should be finally taken to be infinite. The core 
polarization current density becomes 

fcor = -4/ -jj^OiPF ~ \ P \)^-{AV(p)6[V (p)]} (50) 

= " 4 / Wf 6{pF " \P\)f- F {^V(vnVo(p)]} (51) 
= -^/cn^cos^^. (52) 
Then we separate it to several parts as 

fcor=fcor(H)+f cor (F) (53) 
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with 

j c 3 or (#) = -^3 y dn pPF cos9 p 1— 5- (54) 

3cor(F) = -^3 y dn pPF co S p g— (55) 

where AU^ s ^ and AU^^ are shown to be contributions of AC/^) .from Hartree and Fock parts 
of selfenergies, respectively. 

It is not so easy to solve the above equation exactly in the RHF case though it is possible in 
the RH case. On the other hand we have known that the actual momentum dependence is very 
small at least below the Fermi momentum. Then we can suppose that a perturbative way is 
possible with the respect to the momentum dependence. 

Before explaining the actual method, first, we would like to explain the relativistic Hartree 
(RH) case. There the selfenergies are momentum-independent, and the valence current density 
becomes 

Jvar ~nE* F [bb) 

In this case the core-polarization current density is calculated in the following way. 

u z au? 



J cor 



i r 2 n,Ai 

-—3 / dQ p p F cos 6 P — — 

J f 



4 p3p AU? (57) 



3tt 2 E* f 



In the RH calculation 



ml J (2vr)3 "V'E. 



v 



^ f- (58) 



m 2 



Substituting ea.(|58|) into ea.(|57|l. we can get 



3 = 1 PF g 2 4 pp 3 
J 17 E* F X m 2 v 37r 2 E* F iJ 



Jp n o v J^p 

In the RH case the Fermi energy is obtained as 



e F = E F + £l Pb , (60) 



mi 
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and then 



In the low density region below about the saturation, eppd M, so that we can see that the core 
polarization plays a role to cancel the effect of the effective mass in the valence current density 
In the RHF case the contribution from the Hartree part is large, and we cannot use the per- 
turbative way. Since momentum dependence of the selfenergies is not so large, however, the 
difference between Ilo and ITo is small, and then the Hartree part of the total current density 
j 3 (H) can be approximately gotten with the following equation. 

d 3 p ^n 2 (p) 



no(p) 

d 3 p U z (p) 

Mp)- 



(2vr)3 ^IIoG?) 

d 3 P . , r a u z ( P ) 



~^- AU 'I^"°^W> {1 -m h (62) 

where j 3 ar (H) is the valence part of the Hartree current density as 

JvaA H ) ~ 0^-7 — r- (63) 

The space component of the vector selfenergy, which is caused only by the Fock contribution 
in the saturation matter, is very small, and then AU Z is thought to be contributed from the 
Hartree parts as 

AU Z « AC/f = 4/ (64) 

Then the Hartree contribution of the core polarization current density is approximately given 

as 



f cor (H) = fun - f var (H) « l + yH {IS) 3l, ( ri) (05) 

with 



}. (66) 



j(p)" 

As for the Fock part, the momentum dependence of selfenergies are not so large, and its contri- 
bution is not so big in the total current density. Instead of getting it exactly, thus, we can use 
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the perturbative way for the Fock part of the core polarization current density. Along this line 
the variation of the selfenergies are taken to be only the contribution from the valence nucleon 
as 

f S wr^ n{k) Wi^- k) ' <67> 
Af/ »~ -f lw?^wt)^- k) - <68> 

Then we substitute them into the eq. (fS3|) . and get 

jUF) = |4 T3 ^„^«i«^ p - „). (69, 

Next we consider the isovector current. In the similar way we can calculate the isovector current 
density by taking the variation part of the momentum distribution as 

A»(p,r) = ^(p-a). (70) 

In this work the nuclear system is taken to be the isospin symmetric saturated matter plus 
valence nucleon. Thus, the isovector properties can be treated in the perturbative way. Namely 
it can be considered that the Dirac-fields of the valence nucleon isoscalar one, and that those of 
the nucleon in Fermi sea has a very small isovector part coming from the valence nucleon. 
As for the Hartree part we substitute the following Vq (IV) instead of Vq (IS) into the equation 

V H (TV) _ Cl V f d*p 1 Ul(p) 

Vc (IV) -^j T^3-0(P)^{1 - (71) 

As for the Fock part, furthermore, the variations of the selfenergies become 

d 3 k 



. p , , ft f d s k ,,,M*(k) . 

2 J (27ry U (k) 



(p-k), (72) 

*J"M - -rsf J -^3 M*)g^A.(p - k). (73) 
Then the Fock contribution of the isovector core polarization current density is obtained as 

iUF) = f 3 r 3 f ^ Am e p I SMzlSpd±^M AA0 , p _ o) . (74) 

^o(pf) 
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III. RESULTS 

In this section we show results calculated with the above formulation. In this calculation we use 
the parameters (PF1) 



14j for the a- and oj- exchanges to reproduce the saturation properties of 
nuclear matter: the binding energy BE = 16MeV, the incompressibility K = 200MeV and the 
effective mass M*/M = 0.7 at the saturation density po = 0.17fm~ 3 . For comparison we give 
results with momentum-independent selfenergies obtained by the parameter-set PM1 15] that 
gives the same saturation properties. As for the isovector nucleon-nucleon interaction, C% , we 
take the value of PM1. These values are written in Table 1. 

In Fig. ^ we draw the momentum dependence of the scalar selfenergy U s (p) and that of the time 
component of the vector selfenergy Uq(j>). It can be seen that the variation of the momentum- 
dependent selfenergies is only 2.5 % at most below Fermi level, which looks very small. 
In Fig. [21 we show the density-dependence of the Dirac selfenergies U s and Uq on the Fermi- 
surface (a) and the Landau mass (b) with the parameter-sets, PF1 and PM1. Though two 
results of U s and Uq almost agree each other, we can see rather large difference in the Landau 
mass: the value at ps = Po is M^/M = 0.85 in PF1 which is consistent with the value expected 
by the analysis of ISGQR as shown previously. On the contrary, the momentum-independent 
calculation (PM1) gives M£/M = 0.74 which overestimates the excitation energy of ISGQR. 
Hence it is shown that the very small momentum dependence in the nucleon selfenergies enhances 
the Fermi velocity about 15 %, and gives a significant difference in the Landau mass. 
Furthermore we can also see an interesting behaviour of M£ in PF1, namely, its value agrees 
with the bare mass at ps ~ 0.5po an d becomes larger with the decrease of the density. Effects 
of small Dirac effective mass are largely cancelled at low density by the momentum dependence 
created by the one-pion exchange. 

In Fig. |31 we show the density dependence of the isoscalar current density. In the upper column 
(Fig. the solid and chain-dotted lines indicate the total current density and the valence 
current density, respectively. For comparison the current density for the RH approximation 
are also drawn there with the dashed line. From that we can know that the Fock contribution 
suppresses the RH current density, and the core polarization further suppresses it. 
In the lower column (Fig. Eb) we show the contribution from the core polarization. The long 
dashed and dotted lines indicate the core polarization current density contributed from the 
Hartree and Fock parts, respectively. The Hartree contribution reduces the current density, 
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while the Fock contribution enhances it. 

In Fig. 0] we show the isovector current densities; the meaning of each line is the same as that 
in Fig. |3J As for the isovector channel the core polarization does not affect the total current 
density noticeably. 

The most direct observable for the nuclear static current must be the magnetic moment; the 
nuclear medium effect is examined as the discrepancy from the Schmidt value. Thus we should 
compare our results with the normal current density, which is a current density with no medium 
effect and given as 

fo =v F = P -^. (75) 

Ef 

Here we define the following quantity as 

•3 - 3 

Aj 3 = hot ~ Jo _ (?6) 
Jo 

In Fig. [SI we show the density dependence of Aj^. As for the isoscalar current the total curren 
t density, jf ot , almost agrees with the normal current density, jjj. This result is consistent with 
the fact that the core polarization cancels the effect caused by the effective mass. As for the 
isovector, on the other hand, the core polarization does not have significant effect. Here we can 
see very interesting result that the total current density is 10 % less than the normal current 
density in low density region around ps ~ Po/4. This result is consistent with the experimental 
fact that the isovector magnetic moment is 10 % less than the Schmidt value; here we should 
note that the magnetic moment indicates the medium effect in surface region. Of course our 
calculation is performed for the infinite matter, and does not include the contribution from the 
anomalous part. The result does not directly show the experimental observable, but this results 
is very suggestible. 



IV. CONCLUDING REMARKS 



In this paper we have defined a current operator which is consistent with the momentum- 
dependent Dirac fields. It was shown that this current operator automatically include the 
exchange current. This current operator describes the static spatial current which is determined 
by the Landau mass M£ independently of the effective mass M*. This fact is satisfied in all 
cases though we here show it only in the one-pion exchange case which is the most effective. 
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Furthermore we calculate the static current in the system with one valence nucleon on the Fermi 
surface of the saturated nuclear matter. In this calculation we introduce the core polarization 
effect. We can confirm that the core polarization cancels the enhancement caused by the effective 
mass in the isoscalar current, but hardly affects the isovector current. 

As shown in the Ref. jl^ . the very small momentum dependence in the nucleon selfenergies 
enhances the Fermi velocity, even if this momentum dependence is negligibly small for the 
nuclear EOS; in the present calculation the Fermi velocity is enhanced 15 % by the momentum 
dependence caused by the one-pion exchange. Then we succeed to reduce the isovector current 
density in low density region; the value of the current density is almost equivalent to the current 
density without effective mass at ps ~ 0.5/30, and 10% suppressed around ps ~ 0.25po- The 
latter result is consistent with that the observed isovector magnetic moment is 10% smaller than 
the Schmidt value; of course the quantitative conclusion has not been so clear. 
As seen in this paper the momentum-dependent parts, which are non-local in the finite nuclei, 
are very effective in observables related with Fermi velocity even if these parts are small. In 
future we need to discuss effects of the non-local parts of Dirac fields to study nuclear structure 
and reactions. 



The typical value of effective mass is empirically known as M^/M^ = 0.55 — 0.7 [3, li* U, ll* 20]. 
If we use other parameter-set which give smaller effective mass than ours, the effects of the 
momentum-dependent part created by the one-pion exchange does not have a sufficient effect to 
explain the fermi velocity expected from experimental analysis. Exchange forces of other mesons, 
a, u>, r] and 8, also contribute to suppressing the Fermi velocity if we choose the PV-coupling 
for 7r- and rj- nucleon coupling. 

Here we should give a further comment. Bentz et al. have shown in Ref. |ljj] that the Landau 
mass is reduced by the one-pion exchange, which is opposite to ours. This result is consistent 



willi the nonrelativist ie analysis on the magnetic- moment, with the exchange current [^lj- In 
this work Miyazawa showed that the exchange current enhances the convection part, reduces 
the spin part, and totally reduces the isovector magnetic moment. In the calculation [l7| Bentz 
et al. have used the pseudo-scalar (PS) coupling, and the sign of U^ ID was taken to be opposite 
to ours. The full HF calculation with the PS coupling makes too large contribution to the Dirac 
selfenergies [^| while Bentz et al. calculated the Fock term with a perturbative way. Thus a 
calculation with the PV coupling must be more reliable than that with the PS coupling. 
The large discrepancy between the PS and PV coupling comes from relativistic effects in the 
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one-pion exchange. Since the pion mass is smaller than the nucleon fermi energy, relativistic 
effects must be larger in the one-pion exchange. Miyazawa treated the one-pion exchange in the 
nonrelativistic way. Thus it is not strange that our results qualitatively disagree with Miyazawa's 
one. 

In this work we calculate and discuss only the convection current, but not spin current. However 
our formalism is satisfied for the whole Dirac current, thus the spin-current except for the 
anomalous current must be affects by the momentum dependence of the selfenergy in the same 
way. As for the anomalous current we have to consider the vertex correction in another way. 
Since there is no connection between the upper and lower components of the Dirac spinor in 
this current, we can suppose that a perturbative treatment gives sufficient results. It is one of 
the future works. 

11 , , the new current we 



Since the Fock effects are largely seen in the wide energy region 
suggested here also plays an important role with the large momentum-transfer q. In future we 
need to discuss effects of this current in the high momentum transfer phenomena such as the 
quasielastic electron scattering 
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APPENDIX A: ONE-BODY CURRENT OPERATOR 

The density-dependent vertex correction A M is obtained as 

q^{p + q,p) = -Y 1 {p + q)+^{p). (Al) 

Within the one-boson exchange force, the Fock part of the selfenergy is generally written in the 
following way. 

E F (p) = iJ2 C «J ^l a S N (khaA^(p-k) 

-^{{i - K), 7 b ]^(A;)[ 76 , (i> - K)]AW(p - k) (A2) 

where j a ^) is the 7-matrix with the suffix a(b) indicating the scalar, pseudo-scalar, vector, axial- 
vector and tensor, and A^ a ' is the propagator of meson with the quantum number indicated with 
the suffix a. 

Substituting eq. (|A2|) into eq. (jAlj) . we get 

q^(p + q,p) = -£(p + g) + E(p) 

/J At, 
j^- Al a S N {k) la {^ a \p -k + q)- AW(p - k)} 

+ / ^{[{i>-f + ili h ]SN(k)[ lb ,(i>-f + i)]^ b \p-k + q) 

-l(i>-}t), 7 b }S N (k)l lb ,(i>-}t)]A^} (A3) 

When we omit the vertex form factor of the meson-nucleon coupling, the meson propagator is 
given as 

- wh>- < A4 > 

and 

A^(k + q)- A^(k) = -A^ a \k + q)q{2k + q)A i - a \k). (A5) 
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Then the above equation can be rewritten as the following expression, 

/ r Ah, 
j^rSN(kha{A^(p - k + q)(2p -2k + q),A^\p - k)} 

+ <A£c"> / ^{tf,V-l + 4}]s N {k)[v,V-m 

x A^(p -k + q){2p-2k + g) M A^(p - k) 
-[{i>-f + i),l b ]S N (k)h b , l,]^{p -k + q) 

- ky^SNWfay-MAVip-k)} (A6) 
From the above equation we can get the following vertex correction, which we call AW. 

/Ml, 
j^7 a S N (k) la {A^(P -k + q)(2p-2k + q) fl A^(p - k)} 

x A^(p - A; + ?)(2p - 2k + g) M A^(p - fc) 

- - K + rf), 7 6 ]5jv(A;)[7 6 , 7M]A (b) (p -k + q) 

- [ 7 6 ,7M]^(fc)[76,(^-K)]AW(p- fc)} (A7) 

On the other hand the equation (jA2|) can be rewritten as 

E F (p) = iJ2 C «J ^l a S N ( P -kh a A^(k) 

+ J ^[f,l b ]SN{ P -k)[ lb ,f)^ b \k) (A8) 



The eq. JHJ) is given only by the variable transformation (fc — > p — fc) from l|A2|) . 
Substituting eq. (|A8|) into eq. (jAlj) . then, we get 

g^(p + 9 ,p) = + + 

= % Ca j l^^{ S N(p-k + q)-S N (p-k)} 7a A( b \k) 

+ / -^4^7 6 ]{5'iv(p-A ; + g )-S' J v(p-A ; )}[76,K]A( 6 )(A ; ) (A9) 

Using the WT identity (|23|) . the above equation becomes the following expression. 
<fA„ = iq^Caj ^~pS N {k + q)rvS N (k)7 a AW(p-k) 

+ ^£<?6 / ^{l b {i-f)SN{k + q)T,S{k)(i-f) lb A^{p-k) (A10) 
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From the above equation we can get another expression of the vertex correction, which we call 

A( 2 ), as 

A ? } = i J2 C »J i0j-^ aS N{k + q)T,S N {k)la^ a) (p-k) 

+ / ^{[(i>-]t)H b }SN(k + q )r,s(k)[ lb ,(i,-]t)]A( b \ P -k) (Aii) 

The Feynman diagrams corresponding to the approximate identity A^ ~ A^ are shown in 
Fig- El The diagrams in the upper column and in the lower column indicate the vertex corrections 
represented with A*- 2 ) and with A^ 1 ', respectively. This approximation rule exhibits that the 
iteration of the vertex correction is partially equivalent to the meson coupled with the external 
vector field like the meson exchange current. 

The expression of A^ describes the usual vertex correction which is consistent with the Fock 
selfenergy. However its calculation is not easy to be solved because the equation (|A11|) in- 
cludes the iteration scheme. On the other hand we can perform the actual calculation of A^ 
l|A7|) . In addition we should note that the above expression is satisfied in proton and neutron 
independently. Thus we should make the following approximation rule: 

* / ^7 a S N (k + q)T,S N (kh a A^ « 

* / ■^l a S N (k) 7a {A^(p-k + q)(2p-2k + q),A^(p-k)} (A12) 

and 

i J ^{i^-) t + i)^ b ] S NmibAi>-^W b) (p-k + q)(2p-2k + q ),A( b \p-k) 

-[{$-f + i),7 b }S N (k)[lb, 7m]A (6) (p — k + q) 

- [j h , 7 m ]5at(A;)[76, (i> - j£)]A( 6 )(p — k)}. (A13) 

Here we give a comment. If we substitute the full propagator Sn into eq, (|A7|) . we have to 
solve the vacuum polarization, which is also very difficult. In the usual RMF approach we 
usually calculate observables contributed from the nucleon in the Fermi sea by using only the 
density-dependent part Sjj instead of Sp. In the case of the RH, where the selfenergies are 
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momentum-independent, the following equation is satisfied, 
d 4 k 



(2vr)4 
d 4 k 



{S F (k + q)iS D (k)A(p - k) + S D (k)iS F (k - q)A(p -k + q)} 
4 S D (k){A(p-k + q)-A(p-k)}. (A14) 



(27T) 

This equation mentions us that we can describe particle-hole excitations with the usual approxi- 
mation that the density dependent part of the nucleon propagator So © is used in instead 
of the full propagator. The actual momentum dependence is very small, and the equation ()A14|) 
is approximately satisfied in the RHF case, too. From that the rules of eas. (|A12l) and ()A13|) can 
be considered to be available under this approximation. 
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APPENDIX B: ONE-PION EXCHANGE CURRENT OPERATOR 

In this appendix we discuss details of the current operator. As for isoscalar meson-exchanges, 

(2) (2) 

the correction term of the current can be directly used in the above A^ or approximately A^ . 
As for isovector meson exchanges, however, the current operator includes the diagrams of the 
photon connect with the exchange meson (the mesonic current) and that of the photon contact 
with the meson-nucleon vertex (the contact current). 

In the paper we consider only the one-pion exchange, which is considered to be most effective 
because its mass is smallest in meson masses, and the non-locality of the Fock selfenergy with 
the pion exchange is largest. 

The electro-magnetic interaction Lagrangian density is describe as 

C em (x) = Cl m (x) + C™ m {x) + £ c em (x) (Bl) 

with 

Cl m = -eVi(x) 7M i±^V, (B2) 
Cn = -ie{Mx)dnMx)-Mx)d»Mx)}, (B3) 
4 = f {^{x)l^Tiip{x)<j) 2 {x) - ^{x^^^ix^ix) }. (B4) 

Now we separate the vertex corrections to three parts A M = A^ + A m + A£, which is related with 
the above three parts of the electromagnetic interactions. 

Here we assume that the modified Dirac current is contributed only from the proton as 

r M = 7,.^ ZI +A M = f /| l±^ (B5) 
where is an isoscalar operator. The usual vertex correction is given as 

K = J (0 ^ ~ Wr*TaS(k + q)t,^±^S(k)r al5 (i> - J£)A(p - k) (B6) 

Then this equation can be rewritten as 

A; = J,(\ - \t,) (B7) 

with 

^ = ^ f 7?T4 & ~ *HSn(p ~ k)T^S N (k) l5 (i> - ~ k). (B8) 

mi J (2vr) 4 
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Using the rule given by ea. (|A13|) the above equation is approximately written as 
3 fi ~ 3 fi 



ifl f d A k 



fi f 
ill ' 



{#-)£+ &)l$S N {k) lb {$ - ^)A n (p -k + q){2p-2k + q)^A n (p - k) 



- (i> - f + i)-/5SN(kh57^7r(p -k + q) 

- ^S N (k) l5 (i> - K)A^(p - k) } (B9) 

Next we calculate the contribution A m from the so-called pionic current, where the photon 
connects with the pion exchange between nucleons, 

A ™ = " f )(i>-^ + ihBri^-^ + ihBS N (k)^-m—)rMi>-^ 

M J (2vr) 4 ) m n m n 

x (5ii5 2 j - 6 2 iSij)iA(p -k- q)(-i)(2p -2k- q) fl iA(p - k) 

= J (i>-]t + ih 5 s N (kh 5 (i>-]t) 

x A w (p -k + q)(2p-2k + q) ll A 7T (p - k) (BIO) 

and the contribution A c from the co-called contact currents or the Siegel current where the 
photon, pion and nucleon connect at one vertex. These contributions are obtained as 

A m = " / iPk^ C—)^-^ + iHriS N (k)C-^)r j%lti A^(p-k + q) 
^ J (2vr) 4 ) m n m n 

if if 
+( — hiil5rjS N (p)Ti(— )(i> - f)iA^(p -k + q)}(S u 5 2 j - 5 2 i8 lj ) 

= ^fr 3 J ^{ ^-)t + i) l5 S N (kh 5lfl A n (p-k + q) 

+ 757^^(^)75^ - ]i)AAP ~ k) } (Bll) 
From eas. (|B10(l and (|B11|) we can obtain 

A- + A^ = 2r 3 ^. (B12) 

Finally the the vertex correction is given as the summation of the above three contributions, 
which becomes 

A, = A; t + A™ + Kl = |(1 + r 3 )J,. (B13) 



This equation imply that the ansatz shown in eq.(|B5j) is consistent in the present calculation. 
Here we would like to give some comments about the above equation. Under the isospin sym- 
metric matter, first, the above vertex correction satisfies the WT identity as 

l<fj li = Vfr(p + q)-Z$(p) (B14) 
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with 

Z* F (p) = -U[(p)+'fU£(p), (B15) 

where Uf and Uf£ are the Fock contributions of scalar (|19[). and vector (|20|) selfenergies by the 
one-pion-exchange ()19l2(j|) . respectively. 
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9a 9u B a A a 


PF1 


9.699 9.880 27.61 6.134 1.008 20.32 


PM1 


9.408 9.993 23.52 5.651 0.0 20.32 



TABLE I: Parameter sets in this paper. In all cases have used = 138 MeV, m a — 550 MeV, = 783 
McV and C a = 0. 




FIG. 1: Momentum-dependence of the scalar (a) and vector (b) selfenergies. The solid and dashed lines 
indicate the results with PF1 and PM1, respectively. The dotted line denotes the position of the Fermi 
momentum at ps = po- 




FIG. 2: Density-dependence of the Dirac selfenergies U s and Uo on the Fermi-surface (a) and the Landau 
mass (b) the Landau mass (b). The solid and dashed lines indicate the results for PF1 and PM1, 
respectively, and the full square in (b) denotes the value expected empirically from ISGQR. 
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FIG. 3: Density-dependence of the isoscalar nuclear current density (a) and parts of the core polarization 
current density (b). The chain-dotted lines indicate the valence current densities. In the upper pannel (a) 
the dashed and solid lines represent the current density for the RH approximation and the total current 
densities, respectively. In the lower pannel (b) the dashed and dotted line represent the Hartree and Fock 
contributions to the core polarization current densities, respectively. 




FIG. 4: Density-dependence of the isovector nuclear current density (a) and parts of the core polarization 
current density (b). The meanings of lines are the same as those in Fig. EI 
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FIG. 5: The density dependence of the difference between the normal current density and the total current 
density in our model, normalized by the normal current density. The solid and dashed lines represent the 
isoscalar and isovector current densities, respectively. 
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FIG. 6: Fcynman diagram to show the vertex correction for the electro-magnetic current. The upper 
diagrams indicates the usual vertex correction, and the lower ones our approximate current. The solid, 
wave and dashed lines denotes the propagators of nucleon, meson and photon. 



